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Lattice QCD with purely imaginary sources at zero and non-zero temperature. 


1. Introduction 


Exploring a quantum field theory in the presence of external sources is a common way to 
investigate its properties. A relevant example is the response to chemical potentials and background 
fields. We shall consider a class of sources which, in the Euclidean formulation, maintain the 
positivity of the path-integral formulation of the partition function, which for a non-Abelian lattice 
gauge theory is written as 

Z = y = y detM[f/], (1.1) 

so that their effect is directly explorable by numerical Monte-Carlo simulations, without the need 
to deal with a difficult sign problem. 

In fact, such kind of sources are often considered in order to partially circumvent the sign 
problem which is met in the presence of real sources. The best known example is QCD in the 
presence of a baryon chemical potential jJ-s, whose partition function is written as 


Z{T,Hb) = Tr 



' ]■ 


( 1 . 2 ) 


In the path-integral formulation a non-zero Hb amounts to a modification of the quark action den¬ 
sity i//(ni + 7vDv) t/r—^ \j/{m-\-Yv{dy-\-i{Ay — i5ov 1^))^!/ where /r = jUb/ 3 is the quark chemical 
potential. If is real, the anti-hermiticity properties of YvDy are lost and the fermion determinant 
detM(/r) becomes complex (sign problem), so that Monte-Carlo simulations cannot be used, un¬ 
less opposite chemical potentials are given to up and down quarks (isospin chemical potential) or 
to left-handed and right-handed quarks (axial chemical potential /I5). 

Approximate solutions consist in computing Hb derivatives at /ig = 0 (i.e. the cumulants of 
the baryon number distribution, or generalized susceptibilities), in order to reconstruct the free 
energy dependence as a Taylor expansion in jUg ||I[], or in performing importance sampling by 
taking the positive part of the path-integral measure and then reweighting gauge configurations 
with the complex factor |Q]. Both methods have significant limitations: reweighting in general 
fails as the statistics needed to correctly sample the partition function grow exponentially with the 
volume; Taylor expansion is limited to the small /ig region, since the computational effort needed 
to determine the terms of the series grows very rapidly with the order of expansion. 

An alternative is to introduce a purely imaginary quark chemical potential, /i = ijU/ |Q 

[T3| , p3| , in this case the 75-hermiticity of the fermion matrix is preserved 

and the quark determinant is real and positive. In the standard way to introduce a chemical potential 
on the lattice [p^, that amounts to adding a constant 17(1) phase to temporal links appearing in the 
Dirac operator, Ut{n) e‘“^‘Ut{n), i.e. it is equivalent to a twist by an angle dj = Hi/T in the 
fermionic temporal boundary conditions. Charge conjugation implies that the partition function 
is even in /r, so that switching to imaginary is actually like sending and analytic 

continuation can be performed in the T — plane, when possible, to infer properties at real jx 
from simulations at imaginary jX- 

The partition function at imaginary jX contains further interesting information. Z(6/) is of 
course periodic in 6/, with period In, but the additional symmetry of the path-integral measure and 
of the pure gauge term under center transformations, where Nc is the number of colors, implies 
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that the partition function is actually 2n/Nc periodic in 6 /, i.e. 2n periodic in ds = Iin(/iB)/r []T^]. 
Such a periodicity is directly linked to the fact that only globally colorless physical states (hence 
with an integer baryon number ng) contribute to the partition function, both in the low T and in 
the high T phase, however it is differently realized, smoothly or with phase transitions at 0^ = 71 
or odd multiples of it (known as Roberge-Weiss transitions [[1^]), depending on the effective de¬ 
grees of freedom of the theory, which are hadrons below and quarks and gluons above. Lattice 
simulations have provided plenty of evidence of this different behavior from direct simulations at 
imaginary chemical potentials; moreover, the extraction of the Fourier coefficients of the free en¬ 
ergy, considered as a periodic function of 0 ^, gives direct access to the canonical partition function, 
i.e. at fixed baryon number [[ 1 ^. 


A closely analogous example is given by non-Abelian gauge fheories, wifh or wifhouf dynam¬ 
ical fermions, wifh a fopological 0 term. 


Z(0) = 


(1.3) 


where Q = f d'^x-^e^vp^G‘l,yix)G^p„{x) is fhe topological charge, which enters vaiious aspecfs 
of CP-violafion in sfrong inferacfions. Also in fhis case a non-zero 0 makes fhe pafh-infegral 
measure complex, fhus hindering direcf Monle-Carlo simulations. One is usually interested in fhe 
0 -dependence of various quanfifies of physical inferesf, including fhe vacuum energy and, af finite 
T, fhe free energy densify, whose essenfial fealures can be reconsfrucfed, in a Taylor expansion 
approach, in ferms of fhe cumulanfs of fhe probabilify disfribufion of 2 af 0 = 0. The infroducfion 
of a purely imaginary 0 has been considered as an alfernafive since long | 20 ], wifh many interesting 
resulfs obfained in fhe lasf few years for QCD and QCD-like fheories [21, 22, 23, 24, 25]. 

There is an inferesfing analogy befween fhe physics af non-zero, real 0 and fhaf af non-zero, 
imaginary /i, which is linked to fhe 2n periodicify in 0 deriving from fhe facf fhaf Q is integer- 
valued. The analogy involves a sort of dualify: fhe periodicify is non-frivially realized in fhe low 
T phase of fhe fheory, where fhe acfual dependence is on 0 /Nc and phase fransifions appear af 
0 = 71 or odd multiples of if, while a smoofh periodic behavior is realized in fhe deconfined phase. 
However, in fhis case fhe periodic behavior is on fhe same side where a sign problem occurs, so 
fhaf ifs differenf realizations af high T or low T cannof be checked by direcf simulafions, as for 
imaginary /r, buf only indirecfly, fhrough fhe deferminafion of fhe cumulanfs [26]. 


Elecfromagnefic sources have been considered since fhe early days of laffice simulafions fo 
probe fhe elecfromagnefic properties of hadrons [^]. An electromagnetic background field 
modifies the covariant derivative, for a quark with electric charge q, as follows: 


Dp — dpigA‘‘pT‘^ —>■ dp-\- igA^pT‘’ iqap . (1-4) 

The 75 -hermiticity of the Dirac operator, hence the positivity of the quark determinant, is preserved 
only if ap is real valued: on the lattice, SU (3) links appealing in the fermion matrix get multi¬ 
plied by a pure phase factor, Up{n) — Up{n)up{n) where Up{n) ~ exp{iqap{n)) G U (1), which is 
similar to adding an imaginary chemical potential. In the Euclidean path-integral formulation that 
corresponds to a non-zero imaginary electric field, while real electric fields lead to a sign problem. 
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Numerical studies regarding the electric polarizabilities of hadrons usually involve such imaginary 
electric fields: the ground state energy shift, which at the lowest order is proportional to |E 1^ is 
analitically continued from imaginary to real electric fields in order to determine the polarizability. 

Magnetic fields, instead, are real, so that numerical simulations can be performed without any 
technical obstruction. That has opened the way to extensive simulations of QCD in the presence 
of magnetic background fields, whose interesf is justified by the many phenomenological contexts 
where they may play a role, including the physics of compact astrophysical objects like magne- 
tars 1^^, of non-central heavy ion collisions [ ^ ] and of the early Universe [|^], which involve 
magnetic fields going from 10^° Tesla up fo lO'^ Tesla (|eB| ~ 1 GeV^). 

Quantum Field Theories formulated in moving frames, e.g., Lorentz boosted or rotating, can 
also be viewed as systems with external sources, coupled respectively to the total momentum or to 
the angular momentum of the system. For instance, the partition function in a frame moving with 
speed V is written as 


Z(r,v) = Tr (1.5) 

where P is the total momentum of the system. However, the path-integral representation of Eq. ( [Tsl ) 
leads to a sign problem, since the term coupled to v in the Euclidean action is purely imaginary, as 
can be easily verified even in the simplest case of the Euclidean path-integral for a non-relativistic 
particle after a Galileo transformation. Also in this case, the sign problem can be avoided by 
considering the speed v as a purely imaginary quantity, i.e. \ = i^, so that 

Z(r,v) =Tr(e-("-^P)/^) \n). (1.6) 


That, on its turn, can be viewed, for systems which are translationally invariant, as a modification of 
the standard thermodynamical trace, consisting in applying the translation operator exp(/i§ • P)/r to 
the right-hand side eigenstate \n). That has a representation in terms of a path-integral with spatially 
shifted, instead of standard periodic (or anti-periodic), boundary conditions in the time direction 
being the shift). Numerical simulations of such system are interesting, since the -dependence 
of the theory can be used to derive the equation of state [311. 

A similar continuation, from real to imaginary angular velocities, permits the numerical sim¬ 
ulations of QCD in a rotating frame [^, which in this case can be useful and interesting by itself, 
because of the phenomenology related to heavy ion collisions and to compact stars. 


In the following I will cover in more detail only some of the recent developments concerning 
the topics above. In cases where a sign problem exist, many efforts are being dedicated to the 
exploration of new approaches, which could eventually lead to an effective solution. Such efforts 
include Eangevin simulations for generic complex actions [33], lattice simulations on a Eefschetz 
thimble [^, density of states methods [35], formulations in terms of dual variables [36], tensor 
Renormalization Group techniques p^], effective Polyakov loop models [38, 39] and many other 
approaches, which unfortunately will not be discussed in this context. 
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2. Imaginary chemical potentials 


The location of the (pseudo)critical temperature Tc and the nature of the transition as a function 
of jj. are two key issues for the QCD phase diagram. Determining Tc{lJ.) for small values of jj. is 
a well defined goal for analytic continuation: Tc is located for different values of /i/ and then it 
is fitted to a truncated Taylor series or to other expansions (e.g., Fade approximants) and 
continued to real /i. The reliability of the determination of the quadratic term, i.e. the curvature 


Tc 


= l-fC 





( 2 . 1 ) 


has been verified in many differenf ways: by direcf comparison wifh real simulafions for fhe- 
ories wifhouf a sign problem and by comparison wifh fhe oufcome from ofher approaches, 
like reweighfing [^] or Taylor expansion Defermining non-linear confribufions in /r^ is 

harder [^] and faces a fypical problem for analytic confinuafion: confribufions which are sup¬ 
pressed for < 0 may become imporfanf for > 0, because of the different signs in the Taylor 
expansion. 

Recently, such studies have been extended to realistic discretizations of Nf = 2-\-1 QCD with 
physical or almost physical quark masses in two different studies, the first adopting HISQ staggered 
quarks [^, the other adopting stout staggered quarks [15]. In both cases a tree level Symanzik 
improved gauge action and lattices with Nt = 6 and Nt = S have been adopted. Results for the 
curvature fc, Eq. ( |2.1D , are reported in Fig. |] and compared to previous determinations obtained by 
Taylor expansion: analytic continuation brings consistently larger results, with a discrepancy which 
is presently at the 2-a level when systematic effects are taken into account |[I^]. A continuum 
extrapolation of these results is surely needed to clarify the issue, in view of a definite comparison 
with heavy ion phenomenology. 

Like other methods to avoid the sign problem, analytic continuation will hardly say something 
definite regarding the possibile presence of a critical endpoint along the (pseudo)critical line, where 
the transition turns from crossover to first order. Numerical results show that the introduction of 
an imaginary jj. tends to enlarge the region of small quark masses in the Columbia plot where the 
chiral transition is first order 0; then, by analytic continuation (of a critical surface, in this case), a 
real jj. would shrink the same region, thus disfavoring the appearance of a first order at the physical 
point. A safe conclusion, at present, is that the critical point, if any, is not to be found in the small 
jj. region, and may not be directly connected to the chiral first order region taking place at = 0. 
That diminishes the chances to catch it before a complete solution to the sign problem is reached. 

The phase structure at imaginary /r presents many interesting features by itself: the Roberge- 
Weiss first order transition lines, which are present in the high-T region and for Ob = Im(jUB)/r = 
nTi, must stop with critical endpoints at some critical temperature Trw, in order to match with the 
smooth periodic structure at low T. These are the only critical points in the QCD phase diagram 
whose existence can be predicted for sure (even if on the unphysical side); they also coincide, by 
exchange of the temporal direction with a spatial one, with the finite size transition at which charge 
conjugation symmetry breaks © 0- Recent studies, performed both with staggered and 

Wilson fermions [^,10, ^ for A/ = 2 and Nf = 3, have shown that such endpoints turn from first 
to second order, and then to first order again, as the quark masses are tuned from zero to infinity. 
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Figure 1: Curvature of the pseudo-critical line from different studies. From bottom to top: i) and ii) analytic 
continuation, chiral condensate and chiral susceptibility, stout staggered quarks [0; Hi) analytic continua¬ 
tion, chiral susceptibility with /ij = jii, HISQ staggered quarks [|^; iv) Taylor expansion, chiral suscepti¬ 
bility, staggered quark p4-action [p3[]; v) Taylor expansion, chiral condensate, stout staggered quarks [p^. 


with tricritical points separating the different regions. The evolution of this phase structure can 
be monitored from Im(/iB)/r = k down to /Tg = 0, and the use of tricritical scaling can be 
particularly helpful to extrapolate results to the chiral limit. That has permitted to make progress 
on a long-standing issue, like the order of the chiral transition for Nf = 2 QCD in the massless 
limit ID]: the first order region which is present at Im= 7l remains finite down to /Tg = 0, 
in agreement with previous lattice evidence [|4^]. Given the coarse lattices (Nt = 4 with standard 
staggered fermions), such a result is preliminary, but shows that the strategy is promising for future 
improved studies. 

The non-trivial phase structure at imaginary /r has given rise to many model studies and 
to further speculations about possible phenomena related to it [14j. An interesting aspect that we 
would like to stress is the possibility to test lattice methods, currently used to locate the critical 
endpoint (like looking at the radius of convergence of the Taylor series), on the critical points 
present at imaginary /i, whose location is well known from direct simulations; this is something 
that should be done in the future, in order to better assess present systematic uncertainties on the 
determination of the critical point. 


Simulations at imaginary /r provide useful information for many effective models used to 
explore the QCD phase diagram: in the same spirit of analytic continuation, the effective models 
are continued from real to imaginary chemical potentials, where the free parameters of the model 
are fixed against the results of lattice simulations. Examples are given by the effective Polyakov 
line action at finite chemical potential [^], or by studies where a Polyakov-loop-extended Nambu 
Jona-Lasinio (PNJL) model is directly matched to numerical results at imaginary I®. Also 
various derivatives of the free energy with respect to the chemical potentials, i.e. susceptibilities 
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and generalized susceptibilities, can be effectively computed by measuring the response of the 
thermal medium to non-zero imaginary sources, then exploiting analytic continuation. Systematic 
exploratory studies in this direction have been already performed [^. In this context, one should 
also consider alternative ways for inserting the chemical potential on the lattice, see for instance 
Ref. [48], since that could change the effectiveness of analytic continuation. 


3. Electromagnetic background fields 


In the recent past the possibility of performing standard Monte-Carlo simulations, together 
with the interest related to various theoretical and phenomenological issues, has fostered many lat¬ 
tice investigations of the non-perturbative properties of strong interactions in the presence of a mag¬ 
netic background field B. While the initial primary interest, in connection with the phenomenology 
of non-central heavy ion collisions, was in the so-called Chiral Magnetic Effect (CME) [p9|], i.e. 
the electric charge separation taking place along B in the presence of local CP-violating fluctua¬ 
tions, various other topics have been explored, going from magnetic catalysis, i.e. the expected 
enhancement of chiral symmetry breaking induced by B [|5^, to the influence of magnetic fields on 
the QCD phase diagram and the magnetic properties of strongly interacting matter. Some studies 
have provided unexpected results, which have fostered further investigations, both at the level of 
lattice QCD simulations and at that of theoretical models. 

Magnetic catalysis of the QCD vacuum, i.e. the increase of the chiral condensate as a function 
of B, is a general phenomenon which has been confirmed by many investigations, adopting different 
discretizations and number of colors 1011, |5^ , ^ 1^. Instead, less expected developments 


have taken place at finite T. Early determinations [|5^, of the dependence of the pseudo-critical 
temperature Tc on B, obtained by a standard staggered discretization with unphysical quark masses, 
showed a weak increasing behavior of with B, which was also in agreement with most model 
predictions. However, later simulations, obtained with physical quark masses and closer to the 
continuum limit [^], have shown a decreasing behavior of Tc, which goes down by about 20 % 
for cB ~ 1 GeV^; this is typically accompanied, in a region around by a decrease of the chiral 
condensate as a function of B, a phenomenon known as inverse magnetic catalysis. Some evidence 


for inverse magnetic catalysis has been obtained also by simulations of two-color QCD [ |55| ] and by 


adopting dynamical overlap fermions [pSQ. A similar behavior {Tc decreasing with B) was observed 
in the past in the case of chromo-magnetic backgrounds 

Regarding the chiral magnetic effect, lattice QCD studies have focused on various different 
aspects, like: i) looking directly at electric charge separation in a fixed topological and magnetic 
background [Q; ii) looking at the enhancement of longitudinal electric current fluctuations in a 
magnetic background [^] ; Hi) determining the longitudinal electric current in the presence of a 
magnetic background and an axial chemical potential iv) measuring the correlation between 
electric polarization and topological charge in the presence of a magnetic background In 

general, results show a suppression of the effect by about one order of magnitude with respect 
to model predictions (a similar suppression is observed also for the so-called axial magnetic ef¬ 


fect [64]). However, lattice systematics, in particular regarding a correct discretization of chiral 
fermions, should be further checked in the future. 
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Many efforts have been dedicated to investigate the magnetic response of strongly interacting 
matter, i.e. the dependence of its free energy on the magnetic background, which reveals whether 
the material is paramagnetic or diamagnetic. This is, in principle, a clear-cut question for lat¬ 
tice QCD simulations, since it regards the determination of purely equilibrium thermodynamical 
properties, contrary to other properties like the electrical conductivity of strongly interacting mat¬ 
ter 1^5]]. Nevertheless, obtaining a definite answer has been non-trivial since in lattice simulations, 
which usually adopts toroidal (hence compact) geometries in the spatial directions, the magnetic 
field is quantized, so fhaf faking free energy derivatives wifh respecf fo B is nof a well defined 
operation. Affer some firsf affempfs, aimed af defermining only fhe spin componenf of fhe mag- 
nefizafion and suggesfing diamagnetism af all femperafures various improved mefhods have 
been devised. In parficular one can: i) make use of fhermodynamic infegrafion in an exfended 
paramefer space in order fo obfain direcfly free energy differences, insfead of compufing free en- 

j] ; a) compufe fhe magnetization in ferms of pressure anisofropies in fhe 
10 ]; Hi) adopf an exfernal field disfribufion which is nof 


ergy derivatives ||^, | 

direcfion parallel or orfhogonal fo B 
uniform and has zero nef flux across fhe lattice torus, fhus avoiding field quanfizafion and allowing 
for a direcf compufafion of derivafives {Tl\. Consisfenf resulfs have been obfained by fhese mefh- 
ods and by differenl groups, showing fhaf: i) sfrongly inferacfing matter is sfrongly paramagnefic 
above deconfinemenl [ ^ wifh a response which is linear for magnetic fields going 


up to eB ~ 0.1 — 0.2 GeV^ 
in fhe high-T regime [57, 58 


^ and a linear suscepfibilify which grows logarifmically wifh T 
]; ii) fhe magnetic acfivify is instead very weak below Tc', Hi) fhe 
QCD vacuum ifself has a non-linear paramagnefic behavior [70]. A sfrong paramagnefic behavior 
in fhe deconfined phase is also fhe oufcome of many model compufafions [73, [7^ . 

In Fig. § we report a summary of fhe finife T resulfs, which include all sfudies performed 
af or close fo fhe physical poinf (hence excluding early deferminafions wifh Nf = 2 unimproved 
sfaggered quarks [|^], which however provided qualifafively similar resulfs), and are exfrapolafed 
fo fhe confinuum limif in fwo cases. The agreemenf is especially good when resulfs from Refs, [ j^ , 
^ are compared: fhis is expecfed, since all fhree deferminafions adopf fhe same discrefizafion 
of fhe fheory (sfouf sfaggered quarks, while Ref. \12] adopfs HISQ sfaggered quarks), hence fhe 
comparison is a nice fesf fhaf differenf mefhods fo determine fhe magnetic suscepfibilify provide 
consisfenf resulfs. A few issues remain open, regarding in parficular fhe compufafion of non-linear 
corrections fo fhe free-energy dependence of B, which could give imporfanf confribufions fo fhe 
equafion of sfafe in fhe regime of large fields, relevanf to fhe early sfages of fhe Universe, and fhe 
clarification of fhe magnetic properties below T^. Indeed, in a regime where fhe fhermal medium 
can be well approximated by a purely pion gas, compulations based on fhe Hadron Resonance Gas 
(HRG) model [1_4\ lead to fhe predicfion of weak diamagnetic behavior [j^: evidence for such a 
behavior af femperafures around 100 MeV is still marginal [j^, while if clearly emerges when an 


isospin chemical polenfial is added fo fhe syslem, inducing pion condensation 

Lattice simulations have also provided subsfanlial evidence for a direcf influence of magnelic 
fields on fhe gluon sector, mediated by quark loop effecfs, which is expecfed due to fhe sfrongly 
inferacfing nalure of QCD in fhe low energy regime. Large part of magnelic calalysis in fhe QCD 
vacuum is due to dynamical (sea) quark effecls [^, which are also Ihoughl to be af fhe origin of 
inverse magnetic calalysis [76, Magnetic field induced modificafions and anisofropies are also 


visible in observables like fhe gluon acfion densify [p^, or fhe Polyakov loop [76|. Recenfly, 
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Figure 2: Various determinations of the magnetic susceptibility of strongly interacting matter as a function 
of T from studies at or close to the physical point. Results from arXiv: 1309.1142 (HISQ staggered 
quarks) adopt a direct determination of the free energy derivatives in the presence of inhomogeneous mag¬ 
netic fields; results from arXiv:1310.8656 [ |^ (stout staggered fermions) adopt thermodynamic integra¬ 
tion, interpolating in B\ results from arXiv: 1311.2559 [ [^ (stout staggered quarks) have been obtained by 
determining the magnetization via pressure anisotropies; finally, results from arXiv: 1406.0269 (stout 
staggered quarks) adopt thermodynamic integration, interpolating in the quark condensate. 


it has been found that even the static quark-antiquark potential undergoes significat modifications 
in the presence of a magnetic background [77]: the string tension becomes weaker/stronger in 
the direction parallel/orthogonal to B, while the inverse occurs for the Coulomb coupling; such a 
behavior is consistent with some model computations [78]. 

Many studies have approached the problem of determining how the magnetic field influences 
the hadron spectrum. The interest has been driven mostly by the prediction for a superconductive 
phase of strong interactions above a critical magnetic field, which would be signalled by p meson 
condensation 


]. Early numerical evidence about this phenomenon [ [80| ] has not been confirmed 
by later studies [81,|^], however the issue is still open. An open question is also that regarding the 
possible effects of the observed anisotropy of the static quark-antiquark potential on heavy meson 
states, and if these are eventually detectable in heavy ion experiments. 


4. QCD under combined external probes: topological effects in external fields 

In many situations of theoretical or phenomenological interest, one considers two different 
external sources at the same time, in order to probe new aspects of the theory. An illuminating 
example is furnished by the CME itself, where the combination of a magnetic background and 
of an effective local fluctuation of the Q pai'ameter leads to the phenomenon of electric charge 
separation, to be eventually observed in heavy ion experiments. 

Such combinations have been explored by lattice simulations in many contexts. Ref. [ ]6^ has 
explored the combined effect of an axial chemical potential and a magnetic background field in 
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order to obtain a numerical estimate of the CME effect. Ref. |83] has considered a CP-violating 
combination of electric and magnetic backgrounds, i.e. such that E • B 7 ^ 0, in order to determine 
the induced effective Q term in the pure gauge sector, which is linked to the effective QED-QCD 
interactions in the pseudoscalar channel: in this case one relies on a double analytic continuation, 
since both E and the induced effective 0 are purely imaginary. In Ref. [^], an isospin chemical 
potential has been introduced, together with a non-zero magnetic background, in order to explore 
the magnetic properties of isospin unbalanced matter. Various other possibilities could be explored 
along the same line: for instance, were it not for the presence of the sign problem, numerical 
simulation of dense, cold baryon matter in a magnetic background could open the exploration of 
interesting phenomena, like the De Haas - van Alphen effect, which could be relevant to heavy 
astrophysical objects. 

A particularly interesting issue is the study of the combined application of a non-zero 0-term 
and of an electric field. Indeed, a numerical non-perturbative determination of the dependence 
of the electric dipole moment (EDM) of the neutron on 0, at least to the leading linear order, 
dN{^) = 0 + 0 ( 0 ^), would help putting an upper bound on 0 based on the experimental upper 

bound on the EDM of the neutron. This is presently done by estimating by purely dimen¬ 
sional arguments. A possible way to determine r/A?(0) is to measure the three point function of an 
electromagnetic current between two neutron sources at non-zero 0 . (V(pi, 5 i)|V™|V(po,>s'o)) 0 . 
The sign problem emerging from a non-zero 0 can be approached by reweighting configurations 
sampled at 0 = 0 [|^. An interesting alternative is to introduce an imaginary 0 and then perform 
analytic continuation |23], the imaginary 0 being introduced, after an anomalous axial transforma¬ 
tion, as a 75 mass term in the fermionic action; a recent study [25] adopting this strategy has led to 
the upper bound |0| <7.6 x 10^'^. A different strategy, which has been proposed time ago [Q], 
is to introduce an external electric background and measure the EDM from the energy difference 
between states with different spins in the presence of a 0 term: 

(EJ - (EJ = 4')0E, + 0(03). 


Erom a computational point of view, that could be very convenient, since it requires the determi¬ 
nation of a two-point instead of a three-point function. However, one has to face two different sign 
problems at the same time, emerging from the non-zero electric field and from fhe non-zero 0. A 
possibilify could be fo combine imaginary 0 = iBj and imaginary elecfric fields E = /E/, perform 
direcf deferminafions of fhe fwo poinf funcfion 

~ (/)24)0;E;, = -4’0/E/, 


and fhen rely on a double analytic continuation in fhe fwo differenf imaginary sources: preliminary 
resulfs obfained by fhis approach have been reported in Ref. [22]. 


5. Conclusions 

The insertion of external sources or background fields is a common way fo probe fhe properties 
of Quanfum Eield Theories. In fhe confexf of lattice QCD simulations, one needs fhaf fhe posifiv- 
ify of fhe pafh-infegral measure be preserved by fhe insertion, ofherwise sfandard Monle-Carlo 
mefhods cannof be used (sign problem). 
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In the case of magnetic background fields, no such problem appears, so that lattice QCD 
simulations have been able to provide, in the last few years, plenty of information about the QCD 
phase diagram in a magnetic background and the magnetic properties of strongly interacting matter 
at zero and finite T. Peculiar phenomena have also been found, like inverse magnetic catalysis or 
magnetic field induced anisofropies in pure gauge quanfifies, which claim for further invesfigafions 
in fhe nexf few years. 

In ofher cases, which include quark chemical pofenfials, elecfric background fields and fhe 
0-term, fhe inserfion makes fhe pafh-infegral measure complex. A possible way ouf, in fhese cases, 
is fo adopf purely imaginary sources, even if fhaf does nof correspond fo fhe physical inferesfing 
case. One fhen fypically relies on analyfic continuation in order fo exfracf resulfs for fhe real case, 
however, fhe information obfained in fhe presence of imaginary sources can be inferesfing by ifself 
or for direcf comparison wifh ofher models and effective field fheories. 

One of fhe advanfages of performing numerical simulations af non-zero exfernal sources, wifh 
respecf fo fhe sfandard compufafion of expecfafion values af zero source in a Taylor expansion 
approach, is fhaf one fypically needs correlafions funclions and cumulanfs of lower order, wifh 
a greaf benefif in ferms of noise and of compufafional efforf. The frade-off is given by fhe facf 
fhaf one has fo perform dedicafed simulations af non-zero sources, however fhere is still room for 
various improvemenfs fhaf should be explored in fhe fufure. 

In fhis review we have discussed only a small parf of fhe progress made in fhe field in fhe recenf 
pasf. In parficular, we could nof cover fhe large efforts which are being made in order fo perform 
direcf simulafions in fhe presence of sources which make fhe pafh-infegral measure complex, i.e. 
in order fo provide real, insfead of imaginary, solutions fo fhe sign problem. 

References 

[1] C. R. Allton et af, Phys. Rev. D 66, 074507 (2002); Phys. Rev. D 71, 054508 (2005); R. V. Gavai and 
S. Gupta, Phys. Rev. D 68, 034506 (2003); Phys. Rev. D 71, 114014 (2005). 

[2] I. M. Barbour, S. E. Morrison, E. G. Klepfish, J. B. Kogut and M. P. Lombardo, Nucl. Phys. Proc. 
Suppl. 60A, 220 (1998); Z. Eodor, S. D. Katz, Phys. Lett. B 534 (2002) 87; JHEP 0203, 014 (2002). 

[3] M.G. Alford, A. Kapustin, and E. Wilczek, Phys. Rev. D 59, 054502 (1999); M.-P. Lombardo, Nucl. 
Phys. Proc. Suppl. 83, 375 (2000); P. de Eorcrand and O. Philipsen, Nucl. Phys. B 642, 290 (2002); 
Nucl. Phys. B 673, 170 (2003); M. D’Elia and M. P. Lombardo, Phys. Rev. D 67, 014505 (2003); 
Phys. Rev. D 70, 074509 (2004); H. S. Chen and X. Q. Luo, Phys. Rev. D 72, 034504 (2005); 

P. Giudice and A. Papa, Phys. Rev. D 69, 094509 (2004); V. Azcoiti, G. Di Carlo, A. Galante and 
V. Laliena, JHEP 0412, 010 (2004); Nucl. Phys. B 723, 77 (2005); Y. Shinno and H. Yoneyama, 
arXiv:0903.0922 [hep-lat]; L. Karbstein, M. Thies, Phys. Rev. D 75, 025003 (2007); L. K. Wu, 

X. Q. Luo and H. S. Chen, Phys. Rev. D 76, 034505 (2007); P. Cea, L. Cosmai, M. D’Elia and 
A. Papa, JHEP 0702, 066 (2007); S. Conradi and M. D’Elia, Phys. Rev. D 76, 074501 (2007); 

M. D’Elia, E. Di Renzo and M. P. Lombardo, Phys. Rev. D 76, 114509 (2007); K. Nagata and 
A. Nakamura, Phys. Rev. D 83, 114507 (2011). 

[4] P. de Eorcrand and O. Philipsen, JHEP 0701, 077 (2007); JHEP 0811, 012 (2008). 

[5] P. Cea, L. Cosmai, M. D’Elia and A. Papa, Phys. Rev. D 77, 051501 (2008); Phys. Rev. D 81, 094502 
(2010); P. Cea, L. Cosmai, M. D’Elia, C. Manneschi and A. Papa, Phys. Rev. D 80, 034501 (2009); 

P. Cea, L. Cosmai, M. D’Elia, A. Papa and E Sanfilippo, Phys. Rev. D 85, 094512 (2012). 


11 



Lattice QCD with purely imaginary sources at zero and non-zero temperature. 


[6] M. D’Elia and F. Sanfilippo, Phys. Rev. D 80, 014502 (2009); T. Takaishi, P. de Forcrand and 
A. Nakamura, PoS FAT 2009, 198 (2009). 

[7] M. D’Flia and F. Sanfilippo, Phys. Rev. D 80, 111501 (2009); C. Bonati, G. Cossu, M. D’Flia and 
F. Sanfilippo, Phys. Rev. D 83, 054505 (2011). 

[8] P. de Forcrand and O. Philipsen, Phys. Rev. Fett. 105, 152001 (2010). 

[9] O. Philipsen and C. Pinke, Phys. Rev. D 89, 094504 (2014). 

[10] A. Alexandra and A. Fi, PoS FATTICF 2013, 208 (2013). 

[11] F.-K. Wu and X.-F. Meng, Phys. Rev. D 87, 094508 (2013); Phys. Rev. D 90, no. 9, 094506 (2014). 

[12] P. Cea, F. Cosmai and A. Papa, Phys. Rev. D 89, 074512 (2014). 

[13] C. Bonati, P. de Forcrand, M. D’Flia, O. Philipsen and F. Sanfilippo, Phys. Rev. D 90, no. 7, 074030 
(2014). 

[14] K. Nagata, K. Kashiwa, A. Nakamura and S. M. Nishigaki, arXiv: 1410.0783 [hep-lat]. 

[15] C. Bonati, M. D’Flia, M. Mariti, M. Mesiti, F. Negro and F. Sanfilippo, Phys. Rev. D 90, no. 11, 
114025 (2014). 

[16] J. Sugano, J. Takahashi, M. Ishii, H. Kouno and M. Yahiro, Phys. Rev. D 90, no. 3, 037901 (2014); 

J. Takahashi, H. Kouno and M. Yahiro, Phys. Rev. D 91, no. 1, 014501 (2015); J. Takahashi, 

J. Sugano, M. Ishii, H. Kouno and M. Yahiro arXiv: 1410.8279 [hep-lat]; T. Makiyama et al, 
arXiv;1502.06191 [hep-lat]. 

[17] P. Hasenfratz and F. Karsch, Phys. Fett. B 125, 308 (1983); R. V. Gavai, Phys. Rev. D 32, 519 (1985). 

[18] A. Roberge and N. Weiss, Nucl. Phys. B 275, 734 (1986). 

[19] A. Hasenfratz and D. Toussaint, Nucl. Phys. B 371, 539 (1992); S. Kratochvila and P. de Forcrand, 
PoS LAT2005, 167 (2006); Nucl. Phys. Proc. Suppl. 153, 62 (2006); A. Alexandra, M. Faber, 

I. Horvath and K. F. Fiu, Phys. Rev. D 72, 114513 (2005); S. Fjiri, Phys. Rev. D 78, 074507 (2008); 
C. Gattringer and F. Fiptak, Phys. Fett. B 697, 85 (2011); A. Fi, A. Alexandra, K. F. Fiu and 
X. Meng, Phys. Rev. D 82, 054502 (2010); A. Alexandra and U. Wenger, Phys. Rev. D 83, 034502 
(2011); A. Fi, A. Alexandra and K. F. Fiu, Phys. Rev. D 84, 071503 (2011); J. Danzer and 
C. Gattringer, Phys. Rev. D 86, 014502 (2012). 

[20] G. Bhanot and F. David, Nucl. Phys. B 251, 127 (1985). 

[21] V. Azcoiti, G. Di Carlo, A. Galante and V. Faliena, Phys. Rev. Fett. 89, 141601 (2002); M. Imachi, 

M. Kambayashi, Y. Shinno and H. Yoneyama, Prog. Theor. Phys. 116, 181 (2006); B. Alles and 
A. Papa, Phys. Rev. D 77, 056008 (2008); H. Panagopoulos and F. Vicari, JHFP 1111, 119 (2011). 

[22] T. Izubuchi, S. Aoki, K. Hashimoto, Y. Nakamura, T. Sekido and G. Schierholz, PoS FAT 2007, 106 
(2007). 

[23] S. Aoki et al. arXiv;0808.1428 [hep-lat]. 

[24] M. D’Elia and F. Negro, Phys. Rev. Fett. 109, 072001 (2012); Phys. Rev. D 88, no. 3, 034503 (2013). 

[25] F.-K. Guo et al, arXiv: 1502.02295 [hep-lat]. 

[26] C. Bonati, M. D’Elia, H. Panagopoulos andE. Vicari, Phys. Rev. Fett. 110, no. 25, 252003 (2013). 

[27] G. Martinelli, G. Parisi, R. Petronzio and F. Rapuano, Phys. Fett. B 116, 434 (1982); C. W. Bernard, 
T. Draper, K. Olynyk and M. Rushton, Phys. Rev. Fett. 49, 1076 (1982). 


12 



Lattice QCD with purely imaginary sources at zero and non-zero temperature. 


[28] R. C. Duncan and C. Thompson, Astxophys. J. 392, L9 (1992). 

[29] V. Skokov, A. Y. Illarionov and V. Toneev, Int. J. Mod. Phys. A 24, 5925 (2009); V. Voronyuk et al, 
Phys. Rev. C 83, 054911 (2011); A. Bzdak and V. Skokov, Phys. Lett. B 710, 171 (2012); W. -T. Deng 
and X. -G. Huang, Phys. Rev. C 85, 044907 (2012); K. Tuchin, Adv. High Energy Phys. 2013, 490495 
(2013). 

[30] T. Vachaspati, Phys. Lett. B 265, 258 (1991); D. Grasso and H. R. Rubinstein, Phys. Rept. 348, 163 
( 2001 ). 

[31] L. Giusti and H. B. Meyer, JHEP 1111, 087 (2011); JHEP 1111, 087 (2011); L. Giusti and M. Pepe, 
Phys. Rev. Lett. 113, 031601 (2014). 

[32] A. Yamamoto and Y. Hirono, Phys. Rev. Lett. Ill, 081601 (2013) [arXiv:1303.6292 [hep-lat]]. 

[33] E. Karsch and H. W. Wyld, Phys. Rev. Lett. 55, 2242 (1985); G. Aarts and I. O. Stamatescu, JHEP 
0809, 018 (2008); G. Aarts, E. Seiler and I. O. Stamatescu, Phys. Rev. D 81, 054508 (2010); G. Aarts, 
E. Seiler, D. Sexty and I. O. Stamatescu, Phys. Rev. D 90, no. 11,114505 (2014); G. Aarts, 

L. Bongiovanni, E. Seiler, D. Sexty and I. O. Stamatescu, Eur. Phys. J. A 49, 89 (2013); D. Sexty, 
arXiv:1410.8813 [hep-lat]. 

[34] M. Cristoforetti, E. Di Renzo, and L. Scorzato, Phys. Rev. D 86, 074506 (2012); M. Cristoforetti, 

E Di Renzo, A. Mukherjee and L. Scorzato, Phys. Rev. D 88, no. 5, 051501 (2013); M. Cristoforetti e 
al, Phys. Rev. D 89, no. 11, 114505 (2014). 

[35] Z. Eodor, S. D. Katz and C. Schmidt, JHEP 0703, 121 (2007) [hep-lat/0701022]; K. Langfeld, 

B. Lucini and A. Rago, Phys. Rev. Lett. 109, 111601 (2012) [arXiv: 1204.3243 [hep-lat]]; K. Langfeld 
and B. Lucini, Phys. Rev. D 90, no. 9, 094502 (2014) [arXiv: 1404.7187 [hep-lat]]; S. Ejiri, Eur. Phys. 
J. A 49, 86 (2013) [arXiv: 1306.0295 [hep-lat]]; J. Greensite, J. C. Myers and K. Splittorff, JHEP 
1310, 192 (2013) [arXiv: 1308.6712 [hep-lat]]; H. Saito, S. Ejiri, S. Aoki, K. Kanaya, Y. Nakagawa, 

H. Ohno, K. Okuno and T. Umeda, Phys. Rev. D 89, no. 3, 034507 (2014) [arXiv: 1309.2445 
[hep-lat]]. 

[36] C. Gattringer, PoS LATTICE 2013, 002 (2014). S. Chandrasekharan, Phys. Rev. D 82, 025007 (2010); 
J. Bloch, E. Bruckmann, N. Meyer and S. Schierenberg, JHEP 1208, 066 (2012); J. Bloch, Phys. Rev. 
D 86, 074505 (2012) C. Torrero, O. Borisenko, V. Kushnir, B. Alles and A. Papa, PoS LATTICE 
2013, 338 (2014); T. Kloiber and C. Gattringer, arXiv:1410.3216 [hep-lat]; M. Kniely and 

C. Gattringer, arXiv: 1502.00788 [hep-lat]; P. de Eorcrand, J. Langelage, O. Philipsen and W. Unger, 
Phys. Rev. Lett. 113, no. 15, 152002 (2014) [arXiv: 1406.4397 [hep-lat]]; 

[37] A. Denbleyker et ai, Phys. Rev. D 89, no. 1, 016008 (2014); J. Unmuth-Yockey, Y. Meurice, 

J. Osborn and H. Zou, arXiv:141L4213 [hep-lat]. 

[38] J. Greensite and K. Langfeld, Phys. Rev. D 90, no. 1, 014507 (2014); Phys. Rev. D 90, no. 11, 114507 
(2014). 

[39] J. Langelage, M. Neuman and O. Philipsen, JHEP 1409, 131 (2014). 

[40] M. P. Lombardo, PoS LAT 2005, 168 (2006). 

[41] Z. Eodor and S. D. Katz, arXiv:0908.3341 [hep-ph]. 

[42] E. Laermann, E. Meyer and M. P. Lombardo, J. Phys. Conf. Ser. 432, 012016 (2013); E. Laermann, 

E. Meyer and M. P. Lombardo, arXiv: 1304.3247 [hep-lat]. 

[43] O. Kaczmarek, E Karsch, E. Laermann, C. Miao, S. Mukherjee, P. Petreczky, C. Schmidt, 

W. Soeldner and W. Unger, Phys. Rev. D 83, 014504 (2011) [arXiv:1011.3130 [hep-lat]]. 


13 



Lattice QCD with purely imaginary sources at zero and non-zero temperature. 


[44] G. Endrodi, Z. Fodor, S. D. Katz andK. K. Szabo, JHEP 1104, 001 (2011) [arXiv: 1102.1356 
[hep-lat]]. 

[45] T. DeGrand and R. Hoffmann, JHEP 0702, 022 (2007); B. Lucini, A. Patella and C. Pica, Phys. Rev. 
D75, 121701 (2007). 

[46] M. D’Elia, A. Di Giacomo and C. Pica, Phys. Rev. D 72, 114510 (2005); G. Cossu, M. D’Elia, A. Di 
Giacomo and C. Pica, arXiv:0706.4470 [hep-lat]. 

[47] H. Kouno, Y. Sakai, K. Kashiwa and M. Yahiro, J. Phys. G 36, 115010 (2009); Y. Sakai, K. Kashiwa, 
H. Kouno, M. Matsuzaki and M. Yahiro, Phys. Rev. D 79, 096001 (2009); J. Braun, L. M. Haas, 

F. Marhauser and J. M. Pawlowski, Phys. Rev. Lett. 106, 022002 (2011); Y. Sakai, T. Sasaki, 

H. Kouno and M. Yahiro, Phys. Rev. D 82, 076003 (2010); Phys. Rev. D 84, 091901 (2011); 

H. Kouno, M. Kishikawa, T. Sasaki, Y. Sakai and M. Yahiro, Phys. Rev. D 85, 016001 (2012); 

G. Aarts, S. P. Kumar and J. Rafferty, JHEP 1007, 056 (2010); J. Rafferty, JHEP 1109, 087 (2011); 

K. Morita, V. Skokov, B. Friman and K. Redlich, Phys. Rev. D 84, 076009 (2011); K. Kashiwa, 

T. Hell and W. Weise, Phys. Rev. D 84, 056010 (2011); V. Pagura, D. Gomez Dumm and 

N. N. Scoccola, Phys. Lett. B 707, 76 (2012); D. Scheffler, M. Buballa and J. Wambach, Acta Phys. 

Polon. Supp. 5, 971 (2012); K. Kashiwa and R. D. Pisarski, Phys. Rev. D 87, 096009 (2013); 

K. Kashiwa, T. Sasaki, H. Kouno and M. Yahiro, Phys. Rev. D 87, 016015 (2013); T. Misumi and 
T. Kanazawa, JHEP 1406, 181 (2014). 

[48] R. V. Gavai and S. Sharma, Phys. Rev. D 85, 054508 (2012) [arXiv: 1112.5428 [hep-lat]]. 

[49] A. Vilenkin, Phys. Rev. D 22, 3080 (1980); D. Kharzeev and A. Zhitnitsky, Nucl. Phys. A 797, 67 
(2007); D. E. Kharzeev, L. D. McLerran and H. J. Warringa, Nucl. Phys. A 803, 227 (2008); 

K. Fukushima, D. E. Kharzeev and H. J. Warringa, Phys. Rev. D 78, 074033 (2008). 

[50] A. Salam and J. A. Strathdee, Nucl. Phys. B 90, 203 (1975); H. Suganuma and T. Tatsumi, Annals 
Phys. 208, 470 (1991); S. Schramm, B. Muller, A. J. Schramm, Mod. Phys. Lett. A7, 973-982 (1992); 
S. P. Klevansky and R. H. Lemmer, Phys. Rev. D 39, 3478 (1989); K. G. Klimenko, Z. Phys. C 54, 
323 (1992); Theor. Math. Phys. 94, 393 (1993); V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, 
Phys. Rev. Lett. 73, 3499 (1994) [Erratum-ibid. 76, 1005 (1996)]; Nucl. Phys. B 462, 249 (1996); 
Phys. Lett. B 349, 477 (1995). 

[51] P. V. Buividovich, M. N. Chernodub, E. V. Luschevskaya and M. I. Polikarpov, Phys. Lett. B 682, 484 
( 2010 ). 

[52] M. D’Elia and F. Negro, Phys. Rev. D 83, 114028 (2011). 

[53] E.-M. Ilgenfritz, M. Kalinowski, M. Muller-Preussker, B. Petersson and A. Schreiber, Phys. Rev. D 
85, 114504(2012). 

[54] G. S. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. D. Katz and A. Schafer, Phys. Rev. D 86, 071502 

( 2012 ). 

[55] E.-M. Ilgenfritz, M. Muller-Preussker, B. Petersson and A. Schreiber, Phys. Rev. D 89, no. 5, 054512 
(2014). 

[56] M. D’Elia, S. Mukherjee and F. Sanfilippo, Phys. Rev. D 82, 051501 (2010). 

[57] G. S. Bali et al, JHEP 1202, 044 (2012). 

[58] V. G. Bornyakov, P. V. Buividovich, N. Gundy, O. A. Kochetkov and A. Schafer, Phys. Rev. D 90, no. 
3, 034501 (2014). 


14 



Lattice QCD with purely imaginary sources at zero and non-zero temperature. 


[59] P. Cea and L. Cosmai, JHEP 0302, 031 (2003); JHEP 0508, 079 (2005); P Cea, L. Cosmai and 
M. D’Elia, JHEP 0712, 097 (2007). 

[60] M. Abramczyk, T. Blum, G. Petropoulos and R. Zhou, PoS EAT 2009, 181 (2009). 

[61] P. V. Buividovich, M. N. Chernodub, E. V. Luschevskaya and M. I. Polikarpov, Phys. Rev. D 80, 
054503 (2009). 

[62] A. Yamamoto, Phys. Rev. Lett. 107, 031601 (2011). 

[63] G. S. Bali, F. Bruckmann, G. Endrodi, Z. Eodor, S. D. Katz and A. Schafer, JHEP 1404, 129 (2014). 

[64] V. Braguta, M. N. Chernodub, V. A. Goy, K. Landsteiner, A. V. Molochkov and M. I. Polikarpov, 
Phys. Rev. D 89, no. 7, 074510 (2014). 

[65] A. Amato, G. Aarts, C. Allton, P. Giudice, S. Hands and J. I. Skullerud, Phys. Rev. Lett. Ill, no. 17, 
172001 (2013). 

[66] G. S. Bali et al, Phys. Rev. D 86, 094512 (2012) [arXiv; 1209.6015 [hep-lat]]. 

[67] C. Bonati, M. D’Elia, M. Mariti, F. Negro and F. Sanfilippo, Phys. Rev. Lett. Ill, 182001 (2013). 

[68] C. Bonati, M. D’Elia, M. Mariti, F. Negro and F. Sanfilippo, Phys. Rev. D 89, no. 5, 054506 (2014). 

[69] G. S. Bali, F. Bruckmann, G. Endrodi, S. D. Katz and A. Schafer, JHEP 1408, 177 (2014). 

[70] G. S. Bali, F. Bruckmann, G. Endrodi, F. Gruber and A. Schaefer, JHEP 1304, 130 (2013). 

[71] G. S. Bali, F. Bruckmann, G. Endrodi and A. Schafer, Phys. Rev. Lett. 112, 042301 (2014). 

[72] L. Levkova and C. DeTar, Phys. Rev. Lett. 112, no. 1, 012002 (2014). 

[73] M. M. Anber and M. Unsal, JHEP 1412, 107 (2014); V. D. Orlovsky and Y. A. Simonov, Phys. Rev. D 
89, no. 5, 054012 (2014); arXiv: 1406.1056 [hep-ph]; A. N. Tawfik and N. Magdy, Phys. Rev. C 90, 
no. 1, 015204 (2014); K. Kamikado and T. Kanazawa, JHEP 1501, 129 (2015). 

[74] G. Endrodi, JHEP 1304, 023 (2013). 

[75] G. Endrodi, Phys. Rev. D 90, no. 9, 094501 (2014). 

[76] F. Bruckmann, G. Endrodi and T. G. Kovacs, JHEP 1304, 112 (2013). 

[77] C. Bonati, M. D’Elia, M. Mariti, M. Mesiti, F. Negro and F. Sanfilippo, Phys. Rev. D 89, no. 11, 
114502 (2014). 

[78] V. A. Miransky and I. A. Shovkovy, Phys. Rev. D 66, 045006 (2002); E. J. Ferrer, V. de la Incera and 
X. J. Wen, arXiv: 1407.3503 [nucl-th]; R. Rougemont, R. Critelli and J. Noronha, Phys. Rev. D 91, no. 
6, 066001 (2015). 

[79] M. N. Chernodub, Phys. Rev. Lett. 106, 142003 (2011). 

[80] V. V. Braguta, P. V. Buividovich, M. N. Chernodub, A. Y. Kotov and M. I. Polikarpov, Phys. Lett. B 
718, 667 (2012). 

[81] Y. Hidaka and A. Yamamoto, Phys. Rev. D 87, no. 9, 094502 (2013). 

[82] E. V. Luschevskaya and O. V. Larina, Nucl. Phys. B 884, 1 (2014); E. V. Luschevskaya, O. V. Teryaev 
and O. A. Kochetkov, arXiv: 1411.4284 [hep-lat]; E. V. Luschevskaya, O. A. Kochetkov, O. V. Larina 
and O. V. Teryaev, arXiv: 1411.0730 [hep-lat]. 

[83] M. D’Elia, M. Mariti and E Negro, Phys. Rev. Lett. 110, no. 8, 082002 (2013). 

[84] E. Shintani et al, Phys. Rev. D 72, 014504 (2005). 

[85] E. Shintani et al, Phys. Rev. D 75, 034507 (2007). 


15 



